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Recent laboratory experiments have measured fusion cross sections at center- 
of-mass energies low enough for the effects of atomic and molecular electrons 
to be important. To extract the cross section for bare nuclei from these 
data (as required for astrophysical applications), it is necessary to understand 
these screening effects. We study electron screening effects in the low-energy 

X" 

\-4, ' collisions oi Z = 1 nuclei with hydrogen molecules. Our model is based 

on a dynamical evolution of the electron wavefunctions within the TDHF 
scheme, while the motion of the nuclei is treated classically. We find that 
at the currently accessible energies the screening effects depend strongly on 
the molecular orientation. The screening is found to be larger for molecular 
targets than for atomic targets, due to the reflection symmetry in the latter. 
The results agree fairly well with data measured for deuteron collisions on 
molecular deuterium and tritium targets. 
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I. INTRODUCTION AND MOTIVATION 



The determination of nuclear cross sections at the astrophysically most effective energies 
requires in most cases an extrapolation of cross sections that have been measured at the 
lowest possible laboratory energies. It has been generally believed that the uncertainty 
in this extrapolation is reduced by lowering the energies at which data can be taken in 
the laboratory. However, Assenbaum et al. [|| pointed out that this strategy might be 
problematic as, at very low energies, the experimental cross sections do not reflect the 
collisions of bare nuclei. Rather, they are larger due to screening by the electrons present in 
the target and projectile. This screening has to be accounted for in a determination of the 
nuclear cross section required for astrophysical applications. 

Screening by the target electrons has been observed in several low-energy fusion exper- 
iments ( 0] and references therein). Conventionally, the effects of electron screening are 
expressed in terms of an enhancement factor 



where cxbare is the cross section for bare nuclei at center-of-mass energy and ciexp is the 
cross section measured in the presence of target and projectile electrons. The experimen- 
tal analysis generally assumes that the electrons can be thought of as effectively lowering 
the Coulomb energy between the two colliding nuclei by a constant and energy-independent 
energy increment f/g, the screening energy. Furthermore it is assumed that the bare- nuclei 
cross section (usually derived by extrapolation of data at higher energies, which are little af- 
fected by screening) is known. Realizing that Ue « E at those energies currently accessibly 
in the laboratory, the enhancement factor can be written as [||] 



Here, as is conventional in nuclear astrophysics, we have written the cross section in terms 
of the astrophysical S-factor 



fiE) ^exp (E)/abarc(^) 



(1) 
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a(E) = ^exp(-27rr^(i?)), (3) 

where r]{E) = ZiZ2a{fic'^ /2E)^ is the Sommerfeld parameter for nuclei of charges Zi, Z2 and 
reduced mass fi. As is obvious from the exponential energy dependence, electron screening 
effects rapidly become more important with decreasing collision energy. 

The picture underlying the experimental analysis of electron screening effects has been 
confirmed in recent theoretical studies for atomic targets In fact, these calculations 

show that, at those energies at which screening effects have been observed for atomic targets, 
the effective screening potential is nearly constant at distances smaller than the classical 
turning point and can be replaced by a spatially constant screening energy. Furthermore it 
has been found that, to a good approximation, the electrons can be treated adiabatically at 
these low energies, thus allowing one to replace the screening energy by the gain in binding 
energy between the united atom and the asymptotically separated fragment atoms. However, 
it is as yet unexplained why the experimentally determined screening energies apparently 
exceed the adiabatic limits in some cases 0. 

A recent series of impressive experiments has established that the screening energy 

depends on the form of the target (atomic, molecular, solid, etc.). In particular, it has been 
reported that, for the same nuclear reaction, the electron screening is generally smaller for 
molecular targets than for atomic targets. This has been attributed to the fact that, in 
the case of molecular targets, the gain in electron binding is reduced by the energy spent 
to break the molecular bond and by the energy transferred to the spectator nucleus in the 
molecule. Although this dynamical picture qualitatively explains the observed difference 
between molecular and atomic targets, it has not yet been tested theoretically. 

To study screening in molecular targets, we have performed calculations of collisions 
oi Z = 1 nuclei with hydrogenic molecules. In our model the wave functions of the two 
electrons are evolved dynamically within the time dependent Hartree-Fock (TDHF) scheme, 
while the motion of the nuclei is treated classically. That is, we assume that the colliding 
nuclei behave as classical particles beyond the classical turning point, subject to the forces 
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they exert on each other and to the force exerted on them by the electrons and the spectator 
nucleus. However, the electrons are treated quantum mechanically in the TDHF scheme, 
where the wave functions are calculated in the time-dependent potential generated by the 
nuclei and the inner-electron repulsion. In this way, we determine the value of the screening 
potential U{R) between the two coUiding nuclei. Additionally we monitor the force exerted 
on the spectator nucleus by the other charges to test the dynamical picture underlying the 
previous analysis of electron screening effects for molecular targets. 

Our paper is organized as follows. In Section 2, we briefly describe the TDHF method 
applied in our calculation, as well as the numerical methods involved. Our results are 
presented and discussed in Section 3. 

II. THE TDHF METHOD 

We consider the collision of a Z = 1 nucleus with a Z = 1 target nucleus bound in a 
hydrogenic molecule. The nuclei can be treated as point-like and the collision as head-on, 
since the nuclear interaction occurs at distances far smaller than atomic scales. Our choice 
of nuclear coordinates is defined in Fig. 1. Initially, the electron distribution is cylindrically 
symmetric about the axis defined by the two molecular nuclei (z-axis). This symmetry is 
broken as the projectile nucleus approaches. We characterize the trajectory of the projectile 
nucleus by the orientation angle 9 between the initial projectile velocity and the molecular 
axis. In our notation a collision where the projectile passes very near to the spectator nucleus 
is described by 6* 180°. 

With our choice of coordinates the Hamiltonian for the problem reads 

^^=E£+ E J^ + Vp. + Vp. + v.. + v;„ + v„ (4) 

1=1,2 3=P,s,t J 

where m is the electron mass and the Mj are the masses of the projectile (p), target (t) 
and spectator (s) nuclei. Vpt, Vst, and Vps denote the various Coulomb repulsions among the 
nuclei. The Coulomb interaction between electrons (labeled i — 1,2) and the nuclei is 
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and Vee = e^/|ri — is the interaction between the two electrons. 
The Hamiltonian for the isolated hydrogen molecule is given by 

^2 p2 2 

1=1,2^"'- j=s,t'^-^^^3 1=1,2 j=s,t 1^* -^jl 

In the Hartree-Fock approximation to the ground state, the electronic wave function can 
be written as 

vl>(ri,r2,t) = 7A(ri,t)^(r2,t)-^(a(lM2) -/3(l)a(2)) . (7) 

This wave function describes two electrons in a spin-singlet state {a and (3 are the one- 
electron spin states). As our Hamiltonian is independent of the spin coordinates, the spin 
structure of the electron wave function remains unchanged during the collision. 

The initial ground state of the molecule was constructed by evolving a trial wavefunction 
with a fixed nuclear separation R in imaginary time for the molecular Hamiltonian Hmoi , 
which damps out the highest energy components of the wave function. The resulting HF 
ground state is the lowest eigenfunction of the (discretized) Hamiltonian. As long as the 
overlap between the ground state and the trial wavefunction is nonzero, the wave function 
converges to the Hartree-Fock ground state. The nuclear separation has been treated as a 
variational parameter. We find the energy minimum at R = 0.78A, where the converged 
energy is -30.8 eV. Both values are in good agreement with experiment. 

We evolve the spatial wave function ip by solving the TDHF equation 

dip(r, t) 

with 



- $(r,t) + V^ U(r,t) (8) 



$(r t) = + + / ^s^.MM^ (9) 

"^^"'^^ |r-R.(t)r|r-R,(t)r|r-R,(t)| |r - r'l' 

The electron-electron potential has been calculated by solving the appropriate Poisson equa- 
tion subject to a spherical harmonics expansion around the center of electron charge as the 



boundary condition. The TDHF equations have been solved on a three-dimensional grid 
using a second-order expansion in At. The uniform Cartesian grid (with separation 0.26 A) 
had (20,20,40) points in {x, y, z) directions, and if) was assumed to vanish at the boundaries. 

In accordance with our classical treatment of the nuclei, the time dependence of the 
nuclear position vectors Ri(i) is determined from Newton's law, where the force on each 
nucleus is the sum of the Coulomb force of the other nucleus and the force due to the 
electronic charge density. 

III. RESULTS AND DISCUSSION 

The electron response of the target depends on the velocity of the incoming projectile. 
We have studied collisions oia, Z — 1 nucleus with a hydrogenic molecule for various coUision 
angles 9 and for projectile velocities between 0.2q;c and lOac, corresponding to center-of- 
mass bombarding energies E in the d + D2 system between E = 1 keV and 2.5 MeV. In 
our discussion below we will refer to the velocity dependence of the collision in terms of the 
equivalent center-of-mass collision energy E for the d + D2 system; for application to other 
hydrogenic reactions of interest, such as d + T2, these energies must be scaled according to 
the reduced mass of the colliding nuclei. In the calculation the projectile was aimed directly 
at the target nucleus in the molecule, with an initial separation of 6 A. We then followed 
the collision until the projectile had approached the target to the classical turning point. 

In Figs. 2-4 we show "movies" of coUisions at £■ = 2.5 MeV, 25 keV, and 1 keV at 
two selected angles, 9 — 57° and 170°. The time evolution of the electron wave functions is 
represented by contour plots of electron densities, 2\ip{r,t)\'^, projected onto the scattering 
plane spanned by the three nuclei. The positions of the nuclei at the various times were 
indicated by dots. In all cases, the initial electron configuration is the ground state of the 
hydrogen molecule. 

At = 2.5 MeV, the coUiding nuclei move much faster than the average electron velocity. 
In this sudden limit, the projectile falls through the potentials generated by the electrons 
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and the spectator nucleus, adding the respective potential energies to the collision energy. 
As expected, the electron densities at the classical turning point are independent of angle 
(Fig. 5). In the sudden limit, the collision energies are so high that electron screening has 
no influence on the nuclear process. 

In the adiabatic limit, the target and projectile approach each other with a velocity 
much less than the average electron velocity. The electrons are thus able to respond nearly 
instantaneously to the nuclear motion and occupy the energetically most favorable configu- 
ration during the coUision. This situation is illustrated by the coUision at E = 1 keV. Upon 
reaching the classical turning point, the electron configuration corresponds to the ground 
state of a system of charges, Z — 2 and Z — 1, at the appropriate molecular separation. 
Fig. 6 shows the electron configurations at the classical turning point for various collision 
angles. Close inspection shows that the electron configurations vary with angles, so that the 
electrons are not entirely adiabatic even at this low collision energy (see below). 

At = 25 keV (corresponding to a projectile velocity of ac) the nuclear and electronic 
motions are comparable and the electron configuration is strongly distorted as the projectile 
passes through. Furthermore, the electron response depends strongly on the collision angle, 
as can be seen in Fig. 7. 

During the collision we have monitored the expectation value of the forces exerted on 
the three nuclei by the other particles 



\Rpit)-rf \Rpit)-Rtit)\^ \Rpit)-Rtit)\^ ' 

with similar definitions for F{,Fs. The relative force between projectile and target, that is 
exerted by the electrons and the spectator nucleus, along the collision trajectory then is 

F,ei(t) = (F,(t) - F,(t)) . = (F,(t) - F,(t)) . ^ , (11) 

\tXp[t) — iXt[t)\ l-tVrell 

Upon integration of Fj-^i along the collision trajectory we have determined the molecular 
screening potential C/moi, induced by the electrons and the spectator nucleus, 

Umol = / Frel • (iRrel- (12) 



As an example, Fig. 8 shows the relative force and the screening potential for the 
collision aX E = 25 keV and at various angles. Depending on the relative dominance of 
the electrons and the spectator nucleus, the relative force is either attractive or repulsive, 
leading to an oscillatory behavior as a function of target-projectile separation. At large 
angles, the projectile passes through the repulsive Coulomb field of the spectator nucleus on 
its collision trajectory. In these cases, the exerted relative force is repulsive at a few atomic 
units. As expected, the effect is strongest at very large angles. After passing the spectator 
nucleus, the electron cloud dominates; the relative force becomes attractive again. When 
projectile and target are very close, the relative force must vanish. 

The molecular screening potential is attractive in most cases, enhancing the fusion of 
projectile and target. However, at large angles, the screening potential can also become 
repulsive when the projectile passes the spectator nucleus. 

As in 1^ we interpret the screening potential at the classical turning point [Rd) as the 
(negative) screening energy, i.e. {Ue{E,9) = —Ujaoi{Rci,E,9)), which represents the net 
effect of spectator nucleus and electrons on the nuclear fusion process; i.e., we assume that 
the nuclei fuse with an effective energy E + 11^. Fig. 9 shows the screening energy as a 
function of collision energy at various angles. 

The screening energies show a remarkable dependence on the scattering angle. At forward 
angles, Ue{E,6) increases monotonically with decreasing collision energy. However, at back 
angles the screening energy shows a minimum near E = 25 — 30 keV. At extreme back 
angles {6 > 160°) the minimum value even becomes negative, resulting in "anti-screening" 
(or depletion) of the fusion cross section under these kinematical conditions. 

At low energies, the screening potential is attractive at all angles, and for E < 3 keV the 
electronic motions becomes adiabatic and Ue{0) becomes nearly energy independent. How- 
ever, even at these low energies its value depends noticeably on the molecular orientation. 
The screening energy is largest for 9 ^ 90°, while it becomes smallest at large angles. This 
angle-dependence of the screening energy can be understood by considering the effect on 
the colliding nuclei of the electron density between the spectator and the target. At 6 = 90° 
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this electron density induces an attractive force between projectile and target. This force 
becomes less attractive with increasing or decreasing 6 as the force component induced in 
the 2-direction acts like a tidal force, actually hindering the fusing particles. 

In our discussion we have ignored the rotation of the target molecule during the collision. 
This assumption is valid, as the angle AO that the molecule (with rotational energy -Erot) 
would rotate in the time it takes for the projectile (with kinetic energy E) to pass through 
the molecule is small. In particular. 



is only 0.005 radians in the most extreme case discussed here, where £" is 1 keV and E^ot is 



For each screening energy Ue{E,6) one can determine a corresponding enhancement 
factor f{E,6) using Eq. (2). For a comparison with experiment, the enhancement factor 
has to be averaged over angle: 



We have performed this average by 5-point Gaussian quadrature. Again using relation (2), 



Fig. 10. As the extreme backward angles have only little weight in the average, tJe does not 
have a minimum at = 25 — 30 keV as one might expect from Ue{6), but is a monotonically 
decreasing function of energy which changes from the sudden to the (nearly) adiabatic regime 
between 30 keV and 3 keV. 

Fig. 10 also shows the screening energy calculated for a deuteron colliding with an atomic 
D-target. The results are taken from the TDHF calculation of Ref. We observe that the 
screening energy at low collision energies is significantly larger for the D2-target than for 
the atomic D-target. We will argue that this is exceptional. It is the result of the reflection 
symmetry of the d+D system. Asymptotically (large separations) the d+D wave function is 
a 50% mixture of the gerade (positive parity) and ungerade (negative parity) configurations 




(13) 



0.025 eV. 




(14) 



f{E) can be expressed in terms of an average screening energy, Ue{E), which is shown in 
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0. As parity is conserved during the collision, the united atom configuration in the adiabatic 
limit is then 1/2 (He"'"(ls)+He"''(2p)), corresponding to a gain in electronic binding energy 
of 1.5 Ry during the collision. As this symmetry does not hold for the collision with the 
molecular target, the gain is larger and in the adiabatic limit is approximately given by 3 
Ry, corresponding to an electron configuration similar to the one of He"*" (Is). 

More interestingly. Fig. 10 shows that the screening energy is larger in the atomic case 
than in the molecular case for energies E > 20 keV. Obviously this is due to the fact 
that upon angle- averaging the screening in the molecular case is depleted due to the "anti- 
screening" effect when the projectile has to pass the spectator nucleus along its collision 
trajectory. We believe that this result is the main difference between collisions with molecular 
and atomic targets and will hold more generally. One might assume that the collision with 
the molecular target at small 6 roughly imitates the collision on an atomic target. We then 
expect that, in general for colliding systems of non-identical charges the screening energy 
for an atomic target is larger than for a molecular target (see Fig. 9), in agreement with the 
observation reported in [0. We note again that the reffection symmetry in the atomic d+D 
system invalidates our general argument. 

In our calculations we have also examined the force on the spectator nucleus. As an 
example, we consider the collision at E = 25 keV and a.t 6 = 170°. We then calculate a 
gain of kinetic energy of the spectator nucleus of Aitkin = |APp/(4M) ^ 0.1 eV (in the 
lab frame), where AP = JFs{t)dt is the momentum transfer to the spectator. This gain in 
kinetic energy is so small that the spectator does not move before the nuclear fusion process 
happens. We have verified that this is the case for all conditions studied here. Thus, we do 
not support the assumption that the breaking of the molecular bond and the kinetic energy 
transfered to the spectator decreases the screening energy in a collision with a molecular 
target relative to an atomic target 0J^,0. We also find that the assumption that the 
projectile moves on a straight line trajectory is well justified for all orientations appearing 
in the Gaussian quadrature of Eq. (12). 

In Ref. |]lO| the d{d,p) and d{d,n) reactions have been studied down to E = 1.6 keV 
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using a molecular target. Assuming a constant screening energy, the value f/g = 25 ± 5 eV 



has been deduced from the data |T^ between 1.6 keV and ^ 15 keV. As can be seen in Fig. 
10, our calculation does not support the assumption of a constant screening energy in this 
energy regime, but predicts that Ue changes from ^ 20 eV at E = 15 keV to f/g = 40 eV 



at E = 2 keV. Interpreting the screening energy deduced in Ref. [110] as an averaged value 



for the energy interval studied, our calculation appears to be in reasonable agreement with 
data. Fig. 11 compares our prediction for the enhancement of the d + D cross section with 
the data (using Eq. (2) and our calculated Ue{E)), assuming the same parameterized form 
of the bare-nuclei S-factor |jTlI] as in Ref. [|10] . 



Brown et al. [T^ have measured the low-energy t{d, n)a cross section with high precision 



using a T2 target. In Ref. |jT3| it has been argued that the cross section at energies E <1Q 
keV is enhanced due to screening effects. With this assumption and parameterizing the data 
for E = 16 — 70 keV (which are not influenced by electron screening), by a Breit-Wigner 



resonance formula, the authors of [13| were able to consistently describe the Los Alamos 



data |T^ at all energies E <70 keV, a task which previously failed without incorporation of 



screening effects. In Fig. 12 we in fact demonstrate that the t{d, n)a data are well described 
by adding the presently calculated screening energies to the parameterized bare nuclear 
S-factor, as given in ||T3| . 

In conclusion, we have studied the collisions of hydrogen nuclei on hydrogenic molecules 
in a model that evolves the electron distributions dynamically within the TDHF method. 
Our particular interest has been the screening effects of the electrons on the two fusing nuclei. 
We have investigated this effect by determining the molecular screening potential, induced 
by the electrons and the spectator nucleus, on the two colliding nuclei. As in a previous 
study of collisions on atomic targets ||^, we defined the screening energy as the negative 
value of the molecular screening potential at the classical turning point and assumed that 
the screening energy is added to the collision energy of the fusing particle. We found that 
the screening energy depends strongly on the molecular orientation. It is generally smaller 
if the projectile had to pass near the spectator nucleus on its collision trajectory. For the 
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present case of collisions oi Z — 1 nuclei with a hydrogenic molecule we found that the 
screening energy is larger for the molecular target than for the atomic target. This is caused 
by the reflection symmetry in the latter system and we do not expect it to hold in the general 
case. We are planning to investigate this conjecture in future studies of helium collisions on 
hydrogenic targets. 

The work was supported in part by the National Science Foundation, Grant Nos. PHY94- 
12818 and PHY94-20470, and the U.S. DOE at CSUN, Grant No. DE-FG03-87ER40347. 
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FIGURES 

FIG. 1. Coordinates used in our study. The z-axis is defined by the intra molecular separation 
and 9 is the angle between z and the target-projectile separation. 

FIG. 2. Time evolution of the electron density for the collision with a hydrogen nucleus on a 
hydrogenic molecule at a collision energy oi E = 2.5 MeV and at = 57° (left side) and = 170° 
(right side). The electron density is plotted as equally spaced contours. The positions of the nuclei 
are marked by dots. 

FIG. 3. Similar to Fig. 2, but for E = 2h keV. 

FIG. 4. Similar to Fig. 2, but for S = 1 keV. 

FIG. 5. Angular dependence of the electron density when target and projectile are separated 
by a distance corresponding to the classical turning point. The collision energy \s E = 2.5 MeV. 

FIG. 6. Similar to Fig. 5, but for S = 25 keV. 

FIG. 7. Similar to Fig. 5, but ioi E = I keV. 

FIG. 8. Relative force (in eV/A, dashed line) and molecular screening potential (in eV, solid 
line) for a collision of = 25 keV and at selected angles. 

FIG. 9. Screening energy as a function of collision energy and at various angles. 

FIG. 10. Comparison of the angle-averaged molecular screening energy Ue{E) (solid line) with 
the atomic screening energy (dashed line, from Ref. |Q]) as a function of collision energy. 
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FIG. 11. Comparison of the experimental d{d,p)^}i S-factors [10| with the bare-nuclei S-factor 



from Ref. |11| (dashed line) and our estimate in which electron screening is added to the bare- nuclei 



S-factor using Eq. (2). 



FIG. 12. Same as Fig. 11, but for the t{d,n)a reaction. The experimental data are from |12|, 
while the bare-nuclei S-factor is adopted from p^ . 
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